Mode-coupling theory for the slow collective dynamics of fluids adsorbed in 
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We derive a mode-coupling theory for the slow dynamics of fluids confined in disordered porous 
media represented by spherical particles randomly placed in space. Its equations display the usual 
nonlinear structure met in this theoretical framework, except for a linear contribution to the mem- 
ory kernel which adds to the usual quadratic term. The coupling coefficients involve structural 
quantities which are specific of fluids evolving in random environments and have expressions which 
are consistent with those found in related problems. Numerical solutions for two simple models 
with pure hard core interactions lead to the prediction of a variety of glass transition scenarios, 
which are either continuous or discontinuous and include the possibility of higher-order singularities 
and glass-glass transitions. The main features of the dynamics in the two most generic cases are 
reviewed and illustrated with detailed computations. Moreover, a reentry phenomenon is predicted 
in the low fluid-high matrix density regime and is interpreted as the signature of a decorrelation 
mechanism by fluid-fluid collisions competing with the localization effect of the solid matrix. 



I. INTRODUCTION 

Over the last fifteen years, the dynamics of glassform- 
ing liquids under nanoscale confinement has attracted 
much attention. A great variety of molecular systems 
under many different types of confinement has been sub- 
jected to virtually all available experimental techniques, 
while at the same time extensive studies of simple model 
systems using molecular dynamics simulations were un- 
dertaken 

The main reason for this rapidly growing interest is 
that confinement is considered as a potential tool to in- 
vestigate the concept of cooperativity, a key ingredient 
of many glass transition theories [1, [liQ- Indeed, there 
are now many evidences that the dynamics of deeply su- 
percooled liquids is inhomogeneous and that dynamically 
correlated groups of molecules play a crucial role in the 
slowing-down of the dynamics when the temperature is 
decreased. But, up to now, many questions pertaining, 
for instance, to the shape, size or temperature evolu- 
tion of these dynamical heterogeneities remain essentially 
unanswered (for recent progress, see, however, Ref. [^). 

In confinement, geometric constraints associated with 
the pore shape are imposed to the adsorbed fluid and new 
characteristic length scales, like the pore size, come into 
play. Thus, by looking for alterations in the dynamics 
under confinement compared to the bulk, one can hope to 
pinpoint some of the elusive characteristic features of the 
dynamical heterogeneities. For instance, in the simplest 
scenarios, deviations from the bulk behavior are expected 
to occur due to finite size effects, when the typical size of 
the dynamical heterogeneities in the bulk would become 
larger than the pore size. Confinement effects would thus 
provide a ruler to measure dynamical heterogeneities. 

The situation actually turns out more complex. In- 
deed, a direct comparison between the bulk and confined 
fluids is only meaningful if the physical phenomena which 
are specific to confinement have a weak impact on the 
properties of the imbibed fluid or at least if their influence 



is sufficiently well understood that it can be corrected for. 
This is usually not so and strong confinement effects are 
often observed, for instance, the formation of structured 
layers of almost immobile molecules at the fluid-solid in- 
terface. So, in fact, dynamics in conflnement should be 
addressed as a problem of its own, not necessarily with 
reference to the bulk. 

The variety of the systems to consider is immense. 
Porous media can differ in the size, shape and topol- 
ogy of their pore space. They can be made of various 
materials or receive different surface treatments, leading 
to a wide range of fluid-solid interactions which adds to 
the already great variability of the intermolecular inter- 
actions met with usual glassformers. Thus, owing to the 
complexity of the field, a reasonable microscopic theory, 
able to catch at least some of these aspects, could be 
very helpful. Indeed, applied to various models, it would 
allow to explore thoroughly the phenomenology of con- 
fined glassforming liquids and maybe to disentangle the 
contributions of the different physical phenomena which 
interplay in these systems. 

Many porous media, like controlled porous glasses and 
aerogels, are disordered. In the past few years, a very 
useful and quite successful model to deal with this kind 
of systems has been the so-called "quenched-annealed" 
(QA) binary mixture, first introduced by Madden and 
Glandt [§, [lO] . In this model, sketched in Fig. [2 the fluid 
molecules (the annealed component) equilibrate in a ma- 
trix of particles frozen in a disordered conflguration sam- 
pled from a given probability distribution (the quenched 
component). The matrix is assumed to be statistically 
homogeneous, so that, while for any single realization, 
the system lacks translational and rotational invariance, 
all expectation values computed with the matrix prob- 
ability distribution will have the same properties as in 
a truly translationally and rotationally invariant system. 
A common, but not unique, prescription is to take the 
equilibrium distribution of some simple fluid system, so 
that the various samples of the matrix can be thought 
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FIG. 1: Sketch of a QA system. In black, the immobile matrix 
particles. In white, with arrows symbolizing their movement, 
the fluid particles. 



of as the results of instantaneous thermal quenches of 
this original equilibrium system, hence the denomination 
"quenched" for the matrix component. 

Thanks to the assumption of statistical homogeneity, 
as far as the computation of matrix averaged quantities 
is concerned, QA mixtures can be studied with great ease 
using simple extensions of standard liquid state theoreti- 
cal methods. This has been put to good use to derive 
equations describing the structure and the thermody- 
namics of these systems, either via diagrammatic tech- 
niques IqI. Jiol. or by application of the replica trick 

The aim of the present work is now to develop a dy- 
namical theory for QA mixtures, able to deal with the 
problem of the glass transition in confinement. This will 
take the form of an extension to QA systems of the ideal 
mode-coupling theory (MCT) for the liquid-glass transi- 
tion [H, 111, m, [20|. The MCT occupies a central 
place in the study of the dynamics of supercooled liquids 
in the bulk, as one of the very few available microscopic 
theories in this field. It has well known deficiencies, in 
particular the fact that its predicted sharp ergodicity 
breaking transition, the so-called ideal glass transition, 
is always located in the regime of weak to mild super- 
cooling, rather far from the calorimetric glass transition 
point. But, quite remarkably, it has been able to cor- 
rectly predict novel nontrivial relaxation patterns which 
develop in systems like colloidal suspensions with short- 
ranged attractions [2l| or fluids of symmetric dumbbells 
[2^ . when the parameters of these models are varied. It 
seems thus sensible to turn to this theoretical framework 
for a systematic investigation of confined glassforming 
liquids, aiming at understanding, at least qualitatively, 
how the microscopic details of the fluid-solid system im- 
pact its dynamics. Moreover, it is encouraging that re- 
cent computer simulation data for confined fluids could 
be interpreted with the universal predictions of the MCT 
[H, HSj HEI , showing that the mode-coupling glass tran- 
sition scenario might indeed be of some relevance in con- 



finement as well. Note that Ref. 1241 precisely deals with 
a QA system, as do Refs. [26l. l27l. |28|. where other simu- 
lation studies are reported. 

Theories related to the one to be derived have already 
been discussed in the literature. The mode-coupling ap- 
proach to the diffusion-localization transition in the clas- 
sical random Lorentz gas [1^, Is^, HH is of particular rel- 
evance, since this system is actually a QA mixture taken 
in the limit of a vanishing density of the annealed compo- 
nent. Besides, the present theory will borrow some ideas 
from this approach. Other works have dealt with the 
conductor-insulator transition of quantum fluids in ran- 
dom potentials, originally neglecting the fluid-fluid inter- 
actions [3^ . [33I which were later reintroduced by Thakur 
and Neilson [3J|. They appear as special cases of the 
present theory, with additional implicit and uncontrolled 
approximations in the treatment of the effect of random- 
ness on the static fluid correlations [35|. Finally, a sim- 
ilar nonlinear feedback mechanism has been derived for 
the freezing of a polymer chain in a quenched random 
medium using the self-consistent Hartree approximation 
for the Langevin dynamics of the system 36] . 

The paper is organized as follows. In Sec. II, the MCT 
equations for the collective dynamics of a QA mixture are 
derived and discussed, while Sec. Ill explores the relation 
between the dynamics in a self-induced glassy phase and 
in a quenched random environment. In Sec. IV, dynami- 
cal phase diagrams for two simple models are computed, 
which show two basic types of ideal liquid-glass transi- 
tions discussed in more details in Sec. V. Section VI is 
devoted to concluding remarks. Preliminary reports on 
the present work can be found in Refs. [s^l and [ssj . 



II. MODE-COUPLING THEORY FOR 
QUENCHED-ANNEALED SYSTEMS 

In this section, the mode-coupling equations for the 
QA binary mixture are derived and discussed. But, be- 
fore proceeding with the dynamical theory, a few static 
quantities have to be defined. 

As mentioned in the introduction, in a QA system, 
the disordered porous medium is represented by a col- 
lection of Nm rigorously immobile point particles, ran- 
domly placed in a volume V at positions denoted by 
Si,S2, . . . ,sjv^, according to a g iven probability distri- 
bution 7^(81, S2, sat^) [§, Il0l | Its overall density is 
TT-m = Nm/V and the Fourier components of its frozen 
microscopic density, or, in short, its frozen density fluc- 
tuations, are given by 

^™^^g.qs,^ (1) 

where q denotes the wavevector. Their disorder-averaged 
correlation functions, which, because of the assumed sta- 
tistical homogeneity, are diagonal in q and only depend 
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on its modulus q, define the matrix structure factor 



QTiim 



1 



(2) 



where denotes an average over the matrix reahzations. 

The fluid component consists of Nf point particles 
(density nf ~ Nf/V) of mass m, which equilibrate at 
a temperature T in the random potential energy land- 
scape created by the frozen matrix particles. As in the 
bulk, for the present theory, we will be interested in the 
dynamics of the fluid density fluctuations 



Nf 



.7 = 1 



(3) 



where Vj (t) is the position of the fluid particle j at time 
t. At equal times, they allow to deflne the fluid structure 
factor 



1 



and the fluid-matrix structure factor 

1 



(4) 



(5) 



where (• • • ) denotes a thermal average taken for a given 
realization of the matrix, the disorder average being 
performed subsequently. 

For an ordinary binary mixture, the knowledge of the 
above three structure factors would be enough to fully 
characterize the structure at the pair level. This is not 
the case for a QA mixture. Indeed, because the ma- 
trix component is quenched, for any single realization, 
the system lacks translational and rotational invariance. 
It results that, at variance with a bulk fluid, non-zero 
average density fluctuations exist at equilibrium, i.e., 
{Pq) 7^ 0- K is only after averaging over disorder that 

the symmetry is restored, so that (pq) = 0. Thus, one 
is led to consider relaxing and non-relaxing fluid den- 



Non-zero average fluid density fluctuations at equilib- 
rium mean that, even without any dynamical ergodicity 
breaking, they will have time-persistent correlations. In- 
deed, deflning the normalized total density fluctuation 
autocorrelation function 



(p4W/-.(o)) 



(8) 



one expects using standard arguments that 
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This is a general consequence of the fact that the fluid 
evolves in an inhomogeneous environment and we stress 
that this is a true static phenomenon. 

Usually, mode-coupling theories are derived assuming 
that the statics of the problem is solved. For this reason, 
since the calculation of the contribution from the blocked 
correlations is actually a static problem, the theory has 
been developed using the relaxing part of the fluid density 
fluctuations as the central dynamical variable. Attempts 
to derive a dynamical theory starting from the full den- 
sity fluctuations have resulted in complicated equations, 
which appeared unfaithful to the statics of the problem, 
and thus were abandoned. 

The theory is derived using standard projection oper- 
ator methods, as shown in Ref. [11] for bulk systems. As 
usual, in a first step, one obtains a generalized Langevin 
equation for the time evolution of the normalized con- 
nected autocorrelation function of the density fluctua- 
tions 
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(10) 



(bg{t) + n'^cPgit) + n'j dTmq{t-T)4>,{T)^{), (11) 



It is formally the same as for the bulk, i.e. 



sity fluctuations, corresponding to 5p{{t) = pt^{t) - {p{) ^.^^^ ^^^^^^^ conditions 0,(0) = 1, 0,(0) = 0, and 



and (/O^), respectively, and to define the connected fluid 



structure factor 



Sq = j^^{5pk5pL,) 



(6) 



q^keT 
mS'x 



(12) 



and the disconnected or blocked fluid structure factor 

1 



N 



■{pk){pL 



(7) 



such that S^^ ^ Sg + S^. This splitting of the fluid pair 
correlations is well known from the replica theory of QA 
systems, where it leads to the peculiar structure of the 
so-called replica Ornstein-Zernike (OZ) equations [l^.[l3l. 
[3. [isj. which are given in Appendix [K\ for reference. 



The second step involves the calculation of the slow de- 
caying portion of the memory kernel (t) with a mode- 
coupling approach. We will assume that the slow dynam- 
ics is dominated by three types of quadratic variables, 
which can be separated in two classes. In the first class, 
we have variables quadratic in the relaxing density fluc- 
tuations, i5/3k'^Pq_k' i'^ close analogy with bulk MCT p^ . 
In the second class, insp ired by previous studies on the 
Lorentz gas [H, Hflj, |3l[, we consider variables express- 
ing couplings of the relaxing density fluctuations to the 
two frozen density fluctuations present in the problem. 



4 



^PkPq-k ^PkiPci-k)- '^^^ ^^^t variable was omitted 
in all previous works, including a recent account of the 
present theory [sl], with consequences to be discussed 
below. 

The calculation is outlined in Appendix [B] and we only 
quote the result here. It reads mq{t) — rqd{t)+rn'q^'^\t), 
where Tg is a friction coefficient associated with fast dy- 
namical processes, and 



(2^ 



with 



q,k 



(13) 



and 



V 



(1) 

q,k 



q • k ,c q • (q - k) 



•|q-k| 



qC qC 

'-'fc'^|q-k|' 



(14a) 



q-(q-k) 1 

g2 Uf 



qc qb 
'^fc'^|q-k| 



(14b) 

where c^, the Fourier transform of the connected direct 
correlation function, has been introduced. 

Based on these equations, a few general comments 
are in order. First of all, the derived expressions re- 
tain the mathematical structure of the typical mode- 
couplin g eq uations which have been extensively studied 
in Ref. [18l|. Thus, all known properties of the solutions 
of MCT equations, in particular near the transition, ap- 
ply in the case of the QA mixture. A significant addition 
compared to the bulk is however the presence of a linear 
term in the memory function, which opens the possibility 
of continuous ideal glass transitions. 

Like in the theory for the bulk, which is recovered in 
the limit rim — 0, the slow dynamics is fully determined 
by smoothly varying static quantities. Interestingly, no 
explicit reference to the matrix is visible in the equations. 
Indeed, the only required information is n/, S^, Sg, and 
Cg, i.e., quantities characterizing the fluid component of 
the QA system. These functions and the relations be- 
tween them are generically meaningful for the description 
of fluids evolving in statistically homogeneous random 
environments and they are by no means restricted to the 
model of the QA mixture (for the case of non-particle- 
based random fields, see Refs. [H, |4^). Therefore, one 
can expect that the present dynamical theory shares the 
same degree of generality and is applicable in its present 
form out of the strict context of the QA binary mixture. 
In fact, this shows up nicely in the process of deriving the 
equations, since one finds that the contributions resulting 
from the coupling of the random forces to <5pk/3™_k ai'e 
identically zero (see Appendix [B)) . Thus, it is enough to 
consider only two types of quadratic variables, Sp^^Sp^ ^^ 

and <5/3k(Pq-k)' precisely those which are not specific to 
QA systems, to obtain the same dynamical equations. 
Further evidence of the generic character of the present 
equations is given in Appendix [Cl where the dynamics 



of a mean-field spin-glass model in a random magnetic 
field is shown to obey equations with exactly the same 
structure. 

Along the same line, it is noteworthy that the total 
fiuid structure factor S^-^ does not appear in the derived 
dynamical equations. Only 5*^ and 5*^ do. It results 
that, in the framework of the MCT, the global fluid cor- 
relations are of limited relevance to discuss the relation 
between the statics and the dynamics. A much more 
important aspect is rather the balance between the con- 
nected and disconnected contributions to these correla- 
tions. This result is very welcome, since similar differ- 
ences in the roles played by these three types of correla- 
tions are known to be crucial in the physics of QA sys- 
tems. This is best illustrated by the compressibility sum 
rule [ni, dl , which precisely involves and not S^-^ , as 
one would naively expect from the relation for the bulk. 
It is thus reassuring that this feature has been preserved 
despite the uncontrolled approximations involved in the 
derivation of the theory. Moreover, this flnding allows 
one to clarify the issues raised in Refs. [1^, 113] about the 
possibility of a description of the dynamics in confine- 
ment with theories based on structural quantities only. 
Indeed, it was observed there that systems with identical 
global fluid correlations could have significantly different 
dynamics, a result which appeared as a challenge to such 
approaches. The present MCT shows that this is not 
necessarily so and that one should also consider the con- 
nected and blocked correlations before a conclusion can 
be reached. 

About the form of the vertices, one can note the fa- 

(2) 

miliar expression of , which measures the coupling of 

the random forces to Sp(^5p^_^^. It is the same as in a 
bulk system, with connected quantities simply replacing 
the fluid structure factor and direct correlation function. 
This is true for fi^, given by Eq. (|12p. as well. Thus, 
we find that the relaxing part of the density fluctuations 
in a QA system just behaves dynamically like the corre- 
sponding degree of freedom in the bulk. This is not really 
surprising, since a similar correspondence is already visi- 
ble in the statics, for instance in the OZ equation (|A3c[) or 
in the convolution approximation (|B6a[) . which are both 

used in the calculation of . Such a correspondence 
was missing in all our attempts to derive a MCT start- 
ing from the total density fluctuations and this is one of 
the reasons for which the resulting equations were con- 
sidered unsatisfactory. V^^^ is less obvious and combines 



(2) 

features of ^ and of the vertex for the tagged particle 
dynamics in the bulk (Tsj . In particular, like the latter, 
it diverges when q 0. 

Finally, the above equations differ slightly from those 
reported in Ref. [s^l, where a flrst account of the the- 
ory was given. There, 5'|'q_jj| in Eq. (|14b[) is replaced by 



|q-k| ■ 



'"/C|q-k|('^|q-k| 



where is the Fourier trans- 



form of the blocked direct correlation function. The only 
difference between the two approaches is that 5p^{p^ ^^) 
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was not included as a slow variable in the earlier mode- 
coupling scheme. A priori, this is an unjustified approx- 
imation. Indeed, as shown in Appendix [BJ the coupling 
to this variable is actually so strong, that, when both 

'^Pk/°q-k ^^'^ ^PkiPq-k) considered, the contribution 
of the former becomes identically zero, while the effect of 
the random environment is integrally transferred to the 
latter. There are nevertheless circumstances in which the 
absence of Spl_{p{^_^) looks perfectly well motivated, for 
instance when a single particle is moving in the porous 
matrix, corresponding to the limit ^ 0, or when the 
adsorbed fluid is an ideal gas [41)]. In both cases, the 
present MCT reduces to fi^ = g^fc^T/m and 



(2^ 



q-(q-k) 



Mt), (15) 



where is the Fourier transform of the blocked total 
pair correlation function. But, since the only forces ex- 
erted on the fluid are those due to the random matrix, 
on physical grounds, the only mode-coupling contribu- 
tions to the relaxation kernel are expected to come from 
^PkPq-k- previous theories of the Lorentz gas have 
been based on this insight, which leads to the same equa- 
tions as above, except that 'i{'q_k| ^1- fT5|) is replaced 



with h. 



q-k| 



130. 13111. Note that, for an ad- 



sorbed ideal gas, both theories correctly predict that the 
dynamics is independent of n j . 

It seems thus that there is a subtle interplay between 
the approximations involved in the derivation of the 
MCT and the peculiar structure of the static correlations 
in the QA mixture. At present, it is not clear what to 
conclude from this observation, but it is interesting that 
the difference between the two MCT schemes involves 
the blocked direct correlation function Indeed, this 
is a rather delicate object, not easily captured by simple 
approximations. This is best understood in the replica 



framework, where c^{r) is obtained as the zero replica 
limit of the direct correlation function between two non- 
interacting fluid replicas only correlated through their 
common interaction with the matrix [isl lUITsj. It 
results that c^(r) has a highly non-additive character and 
that, for instance, most simple closures of the replica 
OZ equations fail to provide expressions for this func- 
tion which do not vanish identically. Thus, it is a rather 
nontrivial flnding that the MCT approximation scheme 
is actually sensitive to the existence of this function. 



III. SELF-INDUCED CLASSINESS VERSUS 
QUENCHED RANDOM ENVIRONMENT 

From the previous discussion, it is clear that the sepa- 
ration of the fluid density fluctuations into relaxing and 
frozen parts plays a crucial role in the derivation of the 
MCT for QA systems. In that case, the freezing is of 
static origin, due to the random external field generated 
by the quenched matrix. But, as it is well known from 
the MCT for bulk fluids, freezing of the density fluc- 
tuations may also occur dynamically, when the system 
enters in the ideal glassy state. It seems thus interest- 
ing to compare both situations. This can be seen as a 
case of self-induced versus quenched disorder, similar to 
what has been discussed many times in the literature 

The residual relaxation of a bulk fluid in its ideal glassy 
state has been studied in Ref. [47|. It is described by 
mode-coupling equations of the same form as above, ex- 
cept that the characteristic frequency is given by 



and the vertices are 



q Kb 



(16) 
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q • k ^ q • (q - k) 



-Cfc 



-C|q-k| 



{(l-/fe)^a{(l-/|q-k|)^|q-k|} (17a) 



and 



V^'l=nf{{l-U)S,} 



q • k . q • (q - k) , 

— —Ck H 2 ^|q-k| 



{(l-/fe)5,}{/|q-k|5|q_k|} 



(17b) 



where n/ is the density of the fluid, Sq its structure fac- 
tor, Cq the Fourier transform of its direct correlation func- 
tion, and fq the Debye- Waller factor of the glass. Accord- 
ingly, fqSq corresponds to the dynamically frozen part of 
the density fluctuations, while (1 — fq)Sq corresponds to 
their relaxing part. 



The analogies between these equations and Eqs. (fT^ 
and (fT4|) are striking. There is almost a perfect corre- 
spondence between (1 — fq)Sq and 5^ on the one hand, 
and fqSq and on the other hand. Thus, we flnd that, 
irrespective of the mechanism of freezing of the fluid den- 
sity fluctuations, the resulting relaxing and frozen con- 
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tributions essentially play the same role in both systems. 

But there remains one significant difference between 
both sets of equations. Indeed, one of the Fourier trans- 
formed direct correlation functions in Eq. (jl7b[l . the one 
which carries the same wavevector as the frozen part of 
the structure factor, is replaced by a simple constant 
factor in Eq. (fTib]) . The origin of this constant 

term should probably be traced back to the asymmet- 
ric nature of the QA system, where the fluid reacts to 
the matrix but not the other way around, resulting in 
a non-equilibrium character of its static and dynamical 
correlations, in the sense that the matrix is not equili- 
brated with the fluid. It would then reflect the lack of 
dynamical self-consistency between the glassy dynamics 
in a QA system and the frozen background on top of 
which it develops. 

It is tempting to try and obtain the MCT for the QA 
binary mixture as a limiting case of the MCT for the or- 
dinary binary mixture [48l |. where one component would 
become the quenched matrix. The above discussion 
shows that this is not possible. Indeed, based on heuristic 
considerations, one can get close to Eqs. ([TBI) and (fTT]) . 
for instance by canceling all terms which would result in 
a time dependence of the fluid-matrix and matrix-matrix 
correlations, but the non-equilibrium density factor iden- 
tified above seems impossible to generate from the theory 
for the fully annealed system. Thus, as far as the devel- 
opment of theoretical approaches is concerned, and this 
statement is probably not restricted to the MCT frame- 
work, the dynamics of a QA system should definitely be 
considered from the start as different from an infinite 
mass (for Newtonian systems) or a zero bare diffusivity 
(for Brownian dynamics) limit of the dynamics of a fully 
annealed mixture. This is actually not so unexpected 
since, already in the statics, similar difficulties were met 
in early attempts to derive the equations valid for the 
QA mixture starting from those describing fully equili- 
brated systems (see the discussion of Refs. [3, Esf in 
Ref. H). 

Thus, the present theory is not a special case of earlier 
mode-coupling studies of the dynamics of particles mov- 
ing in a glassy matrix 50]. There, the MCT for binary 
mixtures was used in a regime where one component, 
made of big particles, was completely glassy, i.e., both the 
collective and tagged particle correlators did not relax to 
zero, while the other component, made of small spheres, 
was not necessarily localized, i.e., the tagged particle cor- 
relators could relax to zero. In fact, both theories are 
complementary. In the early approach, the glassy ma- 
trix has to be the product of the self-consistent mode- 
coupling dynamics, but the model incorporates thermal 
fluctuations of the solid and aspects of its response to the 
presence of the fluid, while in the present approach, one 
has much more freedom to choose the structure of the 
confining medium, including realistic models of porous 
solids [51| . but the disordered matrix is rigorously in- 
ert. There are anyway qualitative analogies between both 
models. For instance, there is a clear link between the 



fact that the collective motion of the small particles in 
a binary mixture necessarily becomes non ergodic at the 
same point as that of the big particles and the unavoid- 
able existence of static blocked correlations for a fluid in 
a random matrix. 



IV. DYNAMICAL PHASE DIAGRAMS 

We now move to the quantitative predictions of the 
theory, which require numerical solutions of the MCT 
equations. In this section, we report dynamical phase 
diagrams, obtained by mapping, in the parameter space 
of a given model, the domain where lunt^oo (f)q{t) — 0, 
corresponding to the ergodic fluid phase, and the one 
where limt^oo (f'qit) = /? 7^ 0; corresponding to the non- 
ergodic ideal glassy state. The interface between the two 
domains forms the ideal liquid-glass transition manifold; 
in the present work, it will always be a line, since only 
systems with a two-dimensional parameter space are con- 
sidered. Details of the dynamical changes when crossing 
this line are discussed in the next section. 

As already mentioned in Sec. II, a significant difference 
between the MCT equations for QA systems and those 
for bulk glassformers is the presence of a linear term in 
the memory kernel (fT5)l . It results that the discontinuous 
or type B transition scenario known from the bulk, where 
the infinite time limit of (j)q{t) jumps discontinuously from 
zero to a finite value when going from the liquid to the 
glass, is not the only possibility anymore. Continuous or 
type A transition scenarios, where fq grows continuously 
from zero when entering the glassy phase, indeed become 
possible. 

This is best understood by reference to the so-called 
Fin {n > 2) schematic models [l^, in which the time 
evolution of a single correlation function (j){t) is ruled by 
a two parameter memory kernel of the form 

min{t) = vi(p{t) + Vn(l){tT , wi,w„>0. (18) 

It can be readily shown that, when the second term dom- 
inates {vi small), these models have a line of type B tran- 
sitions starting at vi = 0, w„ = n"/ (n— 1)"""'^, and, when 
the first term dominates (u„ small), they have a line of 
type A transitions starting at vi = l,w„ = 0. The study 
of the Fin models also gives us information on the pos- 
sible topologies of the dynamical phase diagrams, which 
depend on the way the transition lines meet. Two simple 
cases are found [l^. If n = 2, the two lines join smoothly 
at a common endpoint, where a topologically stable de- 
generate A3 singularity is located. If n > 2, the two lines 
intersect and, in the glassy domain, only the extension of 
the type B transition line subsists beyond the intersec- 
tion, forming a glass-glass transition line terminated by 
an ordinary A3 singularity. 

It turns out that these two prototypical shapes of phase 
diagrams are obtained with two very simple, closely re- 
lated QA mixture models, to which the present study 
is restricted. In both, the fluid-fluid and fluid-matrix 
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interactions are pure hard core repulsions of the same 
diameter d. The only difference lies in the matrix cor- 
relations. In model I, the matrix configurations are as- 
sumed to be quenched from an equilibrium fluid of hard 
spheres of diameter d, so that the matrix particles do not 
overlap, while in model II the matrix particles are com- 
pletely uncorrelated and overlap freely. In the following, 
both systems will be parametrized by the two dimension- 
less densities 4>f = nufd'^/Q and 0m = Trrim.d^ /6, and the 
Percus-Yevick (PY) approximation [T^ . fisl . [3, 113 ^^^^ 
be used to compute the required structural quantities. 
Note that, in this approximation, c^(r) = 0, so that the 
difflculties mentioned at the end of Sec. II are irrelevant. 

The non-ergodicity parameter fq is a solution of the 
nonlinear set of equations 



1-/. 



(2^ 



(1), 



(19) 



which has to be solved numerically in order to locate the 
liquid and glassy phases when f and 4>m are varied. All 
computations in the present work have been achieved us- 
ing the method of Ref. [1^ , to which the interested reader 
is referred for technical details, and we only provide the 
quantitative information needed to reproduce the present 
results, i.e., that the wavevector integrals have been dis- 
cretized to points on a grid of 128 equally spaced values 
with step size A ~ 0.3835 /d, starting at q^in = A/2 
[S^ l . We have checked by two means that this discretiza- 
tion is not too coarse, in particular with respect to the 
small q divergence of the memory kernel which is cut off. 
First, test calculations at various fluid and matrix densi- 
ties have been performed on a finer q grid, hence with a 
smaller Qmin- Second, in the nf limit, we have com- 
pared our prediction for the diffusion-localization point 
of model I with Leutheusser's analytic result obtained 
within an additional hydrodynamic approximation which 
allows to integrate exactly over the full q range [s^l- In 
all cases, only modest quantitative differences of a few 
percents on the location of the transition points were 
found. 

When dealing with complex transition scenarios, in- 
cluding higher-order singularities and glass-glass transi- 
tion lines, a useful quantity to consider is the largest 
eigenvalue E of the stability matrix of the set of equa- 
tions Indeed, there holds E < 1, and E goes to 
1 when a transition is approached from the strong cou- 
pling side Thus, 1 — E can be used as a convergence 
criterion for the determination of the transition points 
[5^ . In this work, 1 — E < 10~^ was assured for ordinary 
transition points and 1 ~ E < 10~^ was required in the 
regions where the transition lines of types A and B meet. 

The dynamical phase diagrams of models I and II are 
reported in Fig. [21 In the left panel, they are plotted 
in the (0,„,0/) plane. This is the obvious parameter 
space of the problem, but, with this choice of variables, 
no account is given of the differences in structure between 
the two matrix models. This clearly limits the possibil- 
ities of a meaningful comparison between the two sys- 



tems. Thus, in the right panel, both QA mixtures have 
been tentatively parametrized by the same physical con- 
stant, their Henry constant K}j{4'm)- For systems with 
hard core interactions, Knii'm) is equal to the fraction 
of the total volume accessible to the center of an adsor- 
bate particle in a matrix of density 0™ and is probably 
the most simple and generic scalar quantity characteriz- 
ing the confining effect of a solid matrix on an adsorbed 
fiuid. It is given by exp(— 80™) for model II [53|, and by 
exp[— /3/iox(0m)] for model I, where /icx(0m) is the excess 
chemical potential of the equilibrium hard sphere fiuid 
with volume fraction 0™ [58|. The latter has been esti- 
mated followin g th e compressibility route within the PY 
approximation [59|. 

As it can be readily seen in Fig. [21 both phase dia- 
grams essentially have the same overall shape, especially 
when they are plotted as functions of KH{(j)m)- Their 
main qualitative difference at this global scale, which is 
the difference in concavity of their upper parts visible 
in Fig. [2ja), turns out to be rather insignificant, as it is 
the simple and direct consequence of the fact that the 
strength of confinement increases more slowly with 0™ 
in model II than in model I, because of the overlapping 
matrix particles. It completely disappears in Fig. [2fb). 

The nature of the ideal glass transitions met in the 
phase diagrams is just as expected from the study of 
schematic models. Starting from the type B liquid-glass 
transition point for the bulk (at 0„i — 0), where the 
memory kernel is purely quadratic, a line of type B tran- 
sitions develops when increasing 0™, and, starting from 
the type A diffusion- localization point (at 0/ =0), where 
the memory kernel is purely linear, a line of type A tran- 
sitions emerges when increasing 0/. The same is true of 
the way the two lines meet in the phase diagrams. On 
the one hand, for model I, the type A and B lines have 
a common cndpoint, denoted by E, where they smoothly 
join and form a degenerate A3 singularity, like in the 
F12 model. On the other hand, for model II, the two 
lines intersect at a crossing point C and the extension of 
the type B transition line in the glassy domain becomes 
a glass-glass transition line ending with an ordinary A3 
singularity, denoted by E as well. This is the scenario ob- 
tained with the Fin models, for n > 2. Note that, to our 
knowledge, this is the first time that these widely stud- 
ied one equation toy models find physical realizations as 
fluid systems. 

The glass-glass transition line of model II, located be- 
tween points C and E, is barely visible at the scale of 
Fig. [21 It is short and not well separated from the type 
A liquid-glass transition line. Taking into account the 
idealized nature of the predictions of the MCT, it would 
probably be impossible to detect it unambiguously in 
computer simulations of this system. Only the very spe- 
cific features of the dynamics in the vicinity of a higher- 
order singularity, be it degenerate or not, should be vis- 
ible, like lo gari thmic decay laws and subdiffusive behav- 
iors [isl. [2ll.]5^ [60| . But, at least, the present calculation 
shows that this scenario can actually be realized and that 
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FIG. 2: Dynamical phase diagrams of two types of hard sphere quenched- annealed binary mixtures, (a) In the {4>m, 4>f) plane, 
(b) In the {KH{4'm),4'f) plane. Model I: non-overlapping matrix particles. Model II: freely overlapping matrix particles. The 
A3 singularities are indicated with diamonds and denoted by E, while the crossing point between the type A and B transition 
lines in model II is indicated by a cross and denoted by C. 



this does not require any exotic physical ingredient. Now 
that a suitable starting point is available, by playing with 
the parameters of the model, one can try and obtain a 
system for which this glass-glass transition line would be 
extended enough for its signatures to be observable in 
simulations. 

A final requirement for a complete characterization of 
the transitions studied in this work is the knowledge of 
the so-called exponent parameter A [l^. It determines 
many aspects of the dynamics near a transition (see the 
next section) and, for this reason, plays a crucial role in 
the theory. One has < A < 1 and 1/2 < A < 1 for 
type A and B transitions, respectively. Also, A reaches 1, 
its maximum value, at endpoint singularities and jumps 
discontinuously at crossing points. This is thus a useful 
parameter to follow during the computation of the phase 
diagrams. It is plotted in Fig. [3] as a function oi (pf at 
the transition. For model II, points Ci and C2 mark 
the discontinuity associated with the crossing point C 
in Fig. [5] and the existence of glass-glass transitions is 
clearly visible, with values of A given by the line between 
points Ci and E. For the same model, a nonmonotonic 
variation of A along the type B line can be noted as well. 

In addition to these transition scenarios, and formally 
not related to them, another remarkable prediction of 
the present theory is a reentry phenomenon for matrix 
densities higher than the localization threshold (obtained 
for (j)f = 0). Indeed, as shown in Fig. [HJa), for a given, 
not too high (f)m in this domain, ergodicity can be broken 
either by an increase or a decrease of the fluid density. 

In hard core systems, freezing by an increase of the 
fluid density can be qualitatively understood from simple 
free volume arguments, which provide a direct explana- 
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FIG. 3: Exponent parameter A along the transition lines of 
models I and II. The lowest part, where A goes to zero with 
<jf>/, has been omitted for readability. For both models, E 
denotes the A3 singularity, with A = 1. Ci and C2 delimit 
the discontinuity associated with the crossing point between 
the type A and B transition lines in model II. 



tion for the decrease of ^/ at the transition as a function 
of (j),n in the upper part of the phase diagrams. Because 
of the volume excluded by the matrix particles, the larger 
the matrix density is, the smaller the fluid density has to 
be for structural arrest to occur. From the results for 
model I, we might further note that this decrease is such 
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that the total compacity 0/ + (t>m at the transition is a 
decreasing function of (pm as well, reflecting the fact that 
the inclusion of immobile matrix particles in the system 
slows down the dynamics more efficiently than the inclu- 
sion of the same amount of mobile fluid particles. Such 
a behavior, which is hardly surprising, has already been 
observed in molecular dynamics simulations p6l . [27j . 

The possibility of an ergodicity breaking transition by 
a decrease of the fluid density, which is reflected in the 
bottom part of the phase diagrams by the increase with 
0,„ of the transition , is more unexpected. We inter- 
pret this prediction as the signature of a delocalization 
phenomenon induced by fluid-fluid interactions. More 
precisely, we propose that the occasional collisions be- 
tween the fluid particles at low (j)f can destroy the dy- 
namical correlations responsible for the localization of 
individual particles in dense enough matrices. For this 
process to provide an efficient ergodicity restoring relax- 
ation channel, a reasonable criterion is that the localiza- 
tion domains should overlap to allow the fluid particles 
to interact. Thus, the localization length computed at 
(j)f = 0, which decreases when (pm increases, should be 
comparable to the average distance between two fluid 
particles, which decreases when (j)f increases. It is then 
immediate that, starting in the localized state, the larger 
(j)m is, the higher (jjf has to be in order to restore er- 
godicity. The physical implications of this result will be 
discussed in more details in the last section. 

Dynamical scenarios involving reentrant glass transi- 
tion lines, higher-order singularities, and glass-glass tran- 
sition lines, have already been found for colloidal sus- 
pensions with short-ranged attractions [2l| . By analyz- 
ing various contributions to the memory kernel, these 
features were shown to result from the interplay of two 
well defined phenomena, cage effect and bond formation, 
driven by the hard core and attractive parts of the inter- 
action, respectively. It seems thus interesting to attempt 
such an analysis for the present problem. 

Guided by the results of Sec. Ill, we propose that a 
contribution represented by a memory kernel m^'^^°'^ (t) 
derived from mg'^^''(t) by replacing l/uf in Eq. (jl4bp 
with c|q_j,| should be isolated. Indeed, the resulting 
expression then coincides with the one describing the 
residual dynamics of a bulk ideal glass. Therefore, one 
can reasonably expect that, for a QA system, mq'^'^^°\t) 
will provide a fair representation of the mechanism of 
caging by fluid particles in the presence of permanent 
density fluctuations, which is precisely the one at work 
in bulk glassy systems. The remaining linear kernel 
^(conf) j-^^ _ rnf^'~^'^ (t) — ra'q'^^^^ (t) can then be attributed 
to the confinement-specific phenomena, i.e., the local- 
ization effect of the matrix combined with the unusual 
decorrelation mechanism due to fluid-fluid collisions dis- 
cussed above. In principle, one could try to separate 
these two processes, using the fact that, for an ideal 
gas, only the localization effect is present and leads to 
a matrix density at the transition which is independent 



of the fluid density. However, already with the present 
limited separation, the somewhat artificial character of 
the procedure shows up in the form of negative values of 
r7iq™"^^(t), which restrict the density domain where sta- 
ble solutions of the MOT equations can be found, so no 
further decomposition of m'^°^^ (t) was attempted. 

The hypothetical phase diagrams computed with the 
partial kernels mq'^^'^'^\t) and mq'^°'^^\t) for model I are 
reported in Fig. 31 where they are compared to the one 
obtained with the full mg^*^-* (t) . Clearly, the reentry phe- 
nomenon observed in the complete phase diagram can be 
explained by the interplay of the two contributions dis- 
cussed above. But, at variance with the colloidal systems, 
the higher-order singularity is not located in the domain 
where they cross over. The interpretation of this find- 
ing is ambiguous. On the one hand, since a higher-order 
singularity is found on the transition line computed with 
j^^cagc)^^-^ ^ -^yell, one might argue that the singularity 
should be considered as an integral part of the scenario of 
caging by fluid particles. Then, in the relevant domain, 
confinement would simply appear as a modifier of the 
cage structure, progressively changing the nature of the 
ideal jamming transition from discontinuous to continu- 
ous. On the other hand, the bifurcation analysis of the 
MCT scenario shows that a higher-order singularity is 
necessarily formed when two lines of type A and B tran- 
sitions meet. Then, since the type A and B lines arise 
from points representative of systems ruled by confine- 
ment and bulk caging, respectively, one might consider 
that the singularity is the product of the interplay of 
these two phenomena. In favor of the latter interpreta- 
tion, it has been recently suggested that higher-order sin- 
gularities generically result from a competition between 
different arrest mechanisms |6I|]. Note however that the 
conclusions of Ref. (sij are drawn from the consideration 
of bulk systems only, which necessarily enter into the 
glassy state through type B transitions and for which 
there is no mathematical constraint in the theory impos- 
ing a priori the existence of a singularity. It is not imme- 
diate that the proposed statement is valid or even needed 
for systems where type A and B liquid-glass transitions 
coexist. 



V. TYPICAL TRANSITION SCENARIOS 

It results from the previous section that, if one stays 
away from the higher-order singularities, the crossing 
points, and the glass-glass transitions, which require fine 
tuning of the parameters of the models, the present MCT 
for QA mixtures predicts two generic liquid-glass transi- 
tion scenarios. One is discontinuous or type B, the other 
is continuous or type A. In this section, we discuss the 
different features of these transitions which are relevant 
for comparisons of experimental or simulation data with 
the predictions of the theory. The analytic results will 
be quoted without their proofs, which can be found in 
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FIG. 4: Real and hypothetical dynamical phase diagrams of 
model I, computed with the total and partial memory kernels 
rri'^'^\t), rri'q^^''\t) , and m}'q°^^\t). The curves are labeled 
accordingly and the degenerate Az singularities are denoted 
by E. 



Ref. They will be illustrated with detailed compu- 

tations for model I at two matrix densities, = 0.05 
and (pm = 0.22. For the former value, a type B transition 
occurs at </)^ ~ 0.443, with A ~ 0.774. For the latter, a 
type A transition is found at ~ 0.124, with A ~ 0.729. 
For this matrix density, we consider the type A transi- 
tion on the upper branch of the phase diagram, as it is of 
greater relevance for the problem of the glass transition 
in confinement. In both cases, keeping A < 0.8 allows 
one to consider that the transitions are far enough from 
the higher-order singularity. 

Here, a comment on the method of solution of the 
mode-coupling equations is in order. Indeed, as men- 
tioned in the previous section, it involves a cutoff of the 
low q divergence of the memory kernel. Since it has 
been demonstrated that this divergence can change the 
qualitative properties of the solutions in the asymptotic 
regime near the transition [33l [6^ . [63| , this approxima- 
tion and the use of the results of Ref. [l^, which have 
been obtained under the assumption of nonsingular ver- 
tices, might look problematic. In fact, this is not the case 
for physical reasons. Indeed, in the Lorentz gas limit, the 
low q singularity of the memory function has been shown 
to be an ill feature of the mode-coupling approximation 
(6^ . Thus, from a physical point of view, the discretized 
equations considered in the present calculation are actu- 
ally more satisfactory than the continuous ones and the 
results reviewed in Ref. [l3| can legitimately be applied. 

We first consider the density dependence of the non- 
ergodicity parameter /g, shown in Fig. O As it should, 
fq takes a finite value at the ideal glass transition point 
of type B, while it grows continuously from zero at the 



type A transition. In the glassy state, the bifurcation 
analysis of Eq. ([TO]) to leading order yields two universal 
power law behaviors, 

U - fq (X (0/ - d^^f'^ (20) 

for a type B scenario, and 

(X {^f - <j,)) (21) 

near a type A transition. Accordingly, in Fig. [51 the 
curves corresponding to the type B and A transitions 
start with infinite and finite slopes, respectively. 

The wavevector dependence of fq, though not univer- 
sal, is an important prediction of the MCT as well. It 
is reported in Fig. [6l In Fig. [6Ka), fq is shown at the 
type B transition, while in Fig. El^b) , since = at a 
type A transition, the results at <j)f — 1-10/ are plot- 
ted. For reference, the two relevant structure factors Sq 
and Sq at (p'j are also given. Note that in both cases, Sq 
shows maxima both where Sq has maxima or minima. 
At the type B transition, except for the peak at g = 
where fq reaches 1 as a consequence of the diverging ker- 
nel, the non-ergodicity parameter is very similar to the 
one found for a bulk hard sphere fluid and oscillates with 
Sq, which precisely represents bulk-like correlations. The 
overall amplitude is smaller than in the bulk, reflecting 
the fact that, when <f>ra increases, the system evolves to- 
wards a continuous transition scenario. In comparison, 
the non-ergodicity parameter near the type A transition 
appears rather featureless, fq simply decreases from 1 
at g = 0, with, as is increased, a small shoulder de- 
veloping in the wavevector regime where Sq has its main 
peak and Sq its second peak. With the present models 
where the fluid and matrix particles have the same size, 
it is not clear which changes in the static correlations are 
actually responsible for the growth of this contribution. 

Beside these results for the infinite time limit of the 
density correlation functions, the full dynamics is of great 
interest as well. For this computation, which uses the 
algorithm described in Ref. [l^, we follow Ref. [s^ and 
reduce the generalized Langevin equation pip to its form 
valid for Brownian dynamics, 

Tq^q{t)+(j)q{t)+ I dTTJlqit ~ T)^q{T) = Q , (22) 

with 

Tq = U^,,Sy{qdf (23) 

and the initial condition 0^(0) = 1. This simplification 
affects the short time transient part of the dynamics, but 
not its long time properties. In the following, the unit of 
time shall be chosen such that imic = 160. 

The time evolution of the density correlation function 
<j}q{t) aX q c:i 7.09/d, corresponding to the main peak of 
Sq, is reported in Fig. [7] for state points in the vicinity of 
the two transitions discussed above. The curves for other 
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FIG. 5: Fluid density dependence of the non-ergodicity parameter fq = limt^oo <^g(i) for model I at two matrix densities, (a) 
Type B transition for (pm = 0.05. (b) Type A transition for (j>m = 0.22. Data for four representative wavevectors are shown. 
Vertical dotted lines and diamonds indicate the ideal glass transition. 





FIG. 6: Wavevector dependence of the non-ergodicity parameter fq = limt^oo 4'qi't) for model I at two matrix and fluid 
densities, (a) At the type B transition for (fim = 0.05. (b) Near the type A transition for (j)m = 0.22, at (pf — lAtpf- In both 
panels, Sq and Sq at </>^ are reported for reference {Sq has been divided by three for readability). 



values of q are qualitatively similar. In Fig. EJa), the 
two step dynamics typical of the discontimious ideal glass 
transition scenario and well known from the study of bulk 
systems is easily recognized in the curves corresponding 
to the liquid state. The second step, associated to the 
decay from the plateau where (j)q{t) ~ obeys the so- 
called superposition principle, which states that, in this 
regime and for a given q, the shape of the dynamics is 
independent of the state point. The relaxation functions 
only differ through the characteristic time scale Tq , which 



displays a power law divergence. 



r„ cx (0/ - ^'})-\ (24) 



when the transition is approached. One shows that 
where the exponents a and & (0 < a < 1/2, b > 0) are 
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FIG. 7: Time evolution of the connected density correlation function <j)q{'t) at g ~ 7.09/d for model I at two matrix densities, 
(a) In the vicinity of the type B transition for 0™, = 0.05; from left to right, bottom to top: 0/ = 0.99(f)'f, 0.9990^-, 0.999901^, 
O.999990J, 1.OOOO10J, 1.OOO10J, 1.OO10J, l.Ol^J. The horizontal dotted line marks fij, the non-ergodicity parameter at the 
transition, (b) In the vicinity of the type A transition for (;/>„i = 0.22; from left to right, bottom to top: 0/ — 0.9(/)j, 0.99(/)^, 
0.999(/)J, 0.9999^;^, 1.0001</)J, 1.0010}, 1.O10J, l.lcj)}. 



related to A through 



r(i-a)2 _ r(i + 6)2 



r(l-2a) r(l + 26) 



A, 



(26) 



dynamics. For finite values of 
scaling laws 

G{t) = ci3g±{t/T0) 



one finds the 



(29) 



±1 + Xg±{if 



r denoting Euler's gamma function. In the glassy state, 
only the first relaxation step remains and, when t goes 
to infinity, (pqit) reaches fq, which increases with (f)f as 
shown in Fig. ^a,). 

The dynamics near the type A transition visible in 
Fig. EKb) looks significantly different, with a single 
step relaxation scenario, both in the liquid and glassy 
phases. The slowing-down of the dynamics manifests it- 
self through a weak long time tail which extends to longer 

times when 0/ is increased and turns above (j>'j: into a fi- For small t, we have 
nite asymptote which grows as shown in Fig. [5Kb). 

There are nevertheless strong similarities between the 
two dynamics, provided one concentrates on the time do- 
main where 4>q{t) — fq (i-e., where 4>q{t) is small for a 
type A transition). For type B dynamics, this corre- 
sponds to the so-called fast /3 relaxation regime and we 
first specialize the discussion to this case. Then, close 
enough to the transition, a reduction theorem holds, ac- 
cording to which the wavevector and time dependencies 
of (/)g(i) factorize, yielding 



where the master functions g±{t) are the solutions of 



df9±(£-f)5±(f) = 0, (30) 



and the scaling variables obey 



g±(<) = l/£^ 



(32) 



For large t, in the non-ergodic phase, g+{t) goes to a 
constant and Eq. (|20|) results from the expression of C/j. 
In the same time regime, in the ergodic phase, another 
power law behavior sets in for g- (t), yielding the so-called 
von Schweidler decay law. 



5_(£) = -i?? 



(33) 



Mt)^ f^ + Wt)- (27) 
At the critical point, the critical decay law 

Git) = [to/tr (28) 

is obtained, where a is given by Eq. (pS)) and tg is a 
time scale obtained by matching the short and long time 



where b is given by Eq. (|26p and _B is a positive con- 
stant which can be determined by matching Eq. (|33p with 
Eq. (|32p for i; ~ 1. By combining the power law behav- 
iors of Cf3 and in the resulting expression of G(t), one 
recovers the divergence of t^, Eq. ((24)) . 

Moving now to the type A transitions, one finds that, 
both in the ergodic and non-ergodic states, 4>q(t) essen- 
tially behaves near zero as it does when it approaches 
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in the glassy phase in a type B scenario. In particular, 
the critical decay law ([^5]) remains valid at the transi- 
tion. For finite values cj)f — (pj and independently of 
its sign, Eq. (|77|) has to be modified in order to properly 
define G{t) (see Ref. [31 for details), then one finds that 
a scaling law holds, of the form 

G{t) = cg+(t/T), (34) 

where g+{t) is the same function as above. However, 
as in Eq. ([2T|) . the exponents characterizing the scaling 
variables c and r are twice those for a type B transition, 
i.e., 

ccx|0/-0^|, rAo(x |(/./-0jri/^ (35) 

These behaviors are illustrated in Fig. [H where the 
time evolution of \(f>q{t) — f^\ (simply (j)q{t) for the type 
A transition) at g ~ 7.09/d is plotted in a log-log scale in 
order to evidence the power laws and ([55)1 . In this 
graphs, an interesting consequence of the scaling laws 
for small t is clearly visible, which is the symmetric de- 
parture from the critical decay law ([28)1 at long times 
for points located in the liquid and glassy phases at the 
same distance from the transition. This is a particularly 
important feature of the dynamics in the type A scenario. 

The above results do not apply in the vicinity of higher- 
order singularities or crossing points, where a refined 
mathematical analysis is required ^18. , .56. , ,60j . A detailed 
discussion is beyond the scope of the present overview of 
the theory and we shall only mention that such liquid- 
glass transition points signal themselves in the dynam- 
ics through the appearance of logarithmic decay laws 

[illlil,!!!,!!^. 

Finally, already for bulk systems, it is often quite dif- 
ficult to unambiguously demonstrate that the above fea- 
tures are actually present in some experimental or simu- 
lation data, since, as mentioned in the introduction, the 
MCT offers an idealized picture of the glass transition 
phenomenon and there are always significant alterations 
to the theoretical scenario. An additional difficulty can 
be anticipated in the case of QA systems. Indeed, one 
usually has access to the total density correlation func- 
tion (t) and not to the connected function (f)q (t) . Both 
are related through 

= ^Mt) + ^- (36) 

^q 

Thus, the glassy dynamics is modulated by static factors 
and in particular develops itself on top of a state depen- 
dent static background. The separation in the long time 
behavior of 0^(i) of the evolutions which are of purely 
static origin from those which characterize the glassy dy- 
namics might then be delicate. This could be especially 
critical for type A transitions, where the signatures of the 
continuous transition have to be identified just on top of 
the slowly drifting static contribution. Figures demon- 
strating the problem can be found in Ref. ^] . 



VI. DISCUSSION AND CONCLUSION 

In this paper, a mode-coupling theory for the slow 
dynamics of fluids adsorbed in disordered porous solids 
made of spherical particles frozen in random positions has 
been developed. Derived by properly taking into account 
the peculiar structure of the correlations in these systems 
and by including a contribution which had been forgot- 
ten in a previous work [13 , its equations are found to dis- 
play many appealing features. For instance, they show 
universality, in the sense that they do not contain any 
explicit reference to the precise nature of the random en- 
vironment in which the fluid evolves. Also, they compare 
favorably with previous mode-coupling equations derived 
in other contexts, for the residual dynamics in the glassy 
phase of a bulk fluid (Sec. Ill) or for the equilibrium dy- 
namics of a mean-field spin-glass in a random magnetic 
field (Appendix [C)) . Thus, from a formal point of view, 
the theory appears rather satisfactory. 

Nevertheless, a few difficulties remain. First, there is 
the fact that, in the limit of vanishing fiuid-fluid interac- 
tions, the present theory does not coincide with the MCT 
which can be derived by assuming from the start that 
there are no such interactions. Second, there is the di- 
vergence of the memory kernel for small wavevectors and 
the resulting spurious long time anomalies [H, [H, HI] . 
We do not believe that these issues are really harmful, 
even if their handling requires ad hoc approximations, 
but they are worth stressing, since their solutions would 
probably teach us something on the nature of the ap- 
proximations underlying the MCT scheme and on pos- 
sible extensions of the theory. For instance, it has been 
argued by Leutheusser that the inclusion of vertex cor- 
rections within a kinetic theory approach would solve the 
second problem, but no operational scheme was proposed 
for this calculation 62]. 

The numerical solution of the MCT equations for two 
simple fluid-matrix models leads to a variety of transition 
scenarios, which are either discontinuous for dilute matri- 
ces or continuous for dense matrices. Depending on the 
model, in the intermediate region where the nature of the 
transition changes, degenerate or genuine higher-order 
singularities and glass-glass transition lines are found. 

Another remarkable prediction of the theory is the pos- 
sibility of a reentry phenomenon for high matrix densities 
above the localization threshold, which has been inter- 
preted as the signature of a decorrelation process induced 
by fluid-fluid collisions. 

Before going further, one should note that, strictly 
speaking, this prediction of the theory cannot be cor- 
rect in the case of hard core fluid-matrix interactions. 
Indeed, as recently confirmed_by extensive computer sim- 
ulations of the Lorentz gas [6a ], the localization transi- 
tion is driven by the percolation transition of the ma- 
trix void space, i.e., localization occurs because, above a 
certain critical matrix density (j)'^^ the void space only 
consists of finite disconnected domains. In such a sce- 
nario, it is obvious that, whatever the fluid density, it 
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FIG. 8: Critical dynamics of the connected density correlation function 4iq{t) at g ~ 7.09/d for model I at two matrix densities, 
(a) In the vicinity of the type B transition for (jtm ~ 0.05. As it should, only the data for (j)f < 4>'} appear in the bottom 
panel (note its reverted y scale, in order to make it graphically clear that, despite the use of the absolute value, a decay law 
is represented), (b) In the vicinity of the type A transition for <j)m = 0.22. See Fig. [7] for attribution of the curves. The 
analytically derived power law behaviors (|32p and (|33p are shown as dashed straight lines. 



is impossible to have an ergodic system above c/)^, since 
fluid-fluid interactions will never change the geometry of 
the matrix. This contradiction between the percolation 
theory and the MCT clearly raises the issue of the rela- 
tion between the two approaches, for which we propose 
the following simple argument. The MCT applied to the 
problem of the diffusion-localization transition attempts 
to capture the onset of percolation in an indirect way, by 
giving an account of the increasingly correlated nature of 
the fluid-matrix collision events as the threshold is ap- 
proached. In this respect, one should note that none of 
the static structure functions on which it is based does 
show a sensitivity to the phenomenon of percolation. Fol- 
lowing Leutheusser, the theory works at the level of a self- 
consistent treatment of ring collision processes j62| . This 
turns out to be enough to predict a diffusion-localization 
transition, but, clearly, the infinite sequences of corre- 
lated collisions which would really reflect the permanent 
trapping of the fluid particle in a flnite domain above 
the percolation threshold are missing. From this incom- 
plete characterization of the dynamical processes escort- 
ing the percolation phenomenon, it results that the MCT 
diffusion-localization transition is actually an ideal ver- 
sion of the true percolation transition, in the usual sense 
that the MCT predicts ideal glass transitions, and that 
the theory is not able to detect that, in an exact treate- 
ment, the percolation threshold flxes an absolute limit 
to diffusive behavior. When fluid-fluid collisions come 
into play at flnite fluid densities, this leaves room for the 
prediction of a reentry phenomenon in contradiction with 
percolation theory. We believe that it is for the same rea- 
son that the MCT also misses the fact that, at any matrix 



density, there is always a non-vanishing probability that 
particles will be trapped in a flnite domain disconnected 
from the rest of the void space, so that the exploration of 
the available void space is never completely ergodic [6^ . 

It is thus clear that the prediction of a reentrant be- 
havior of the ergodicity breaking transition line in the 
low fluid density domain should not be taken too lit- 
erally. In fact, a reasonable expectation based on this 
flnding is that, below the localization threshold, but in 
the regime where transient trapping effects are impor- 
tant, there might be an acceleration of the dynamics due 
to fluid-fluid collisions. Interestingly, such a behavior has 
already been observed in a computer simulation study of 
a two dimensional lattice gas model with flxed randomly 
placed scatterers [67*1 . Indeed, it was found that, starting 
from the zero fluid density limit, the diffusion coefficient 
of a tagged particle flrst increases with the fluid density. 
An interpretation in terms of a decorrclation process sim- 
ilar to the one discussed in Sec. IV was then proposed 
and validated by varying the parameters of the model. 

As a possible source of the difficulties of the theory, 
one might blame the fact that it works at the level of dis- 
order averaged quantities. Indeed, the procedure of aver- 
aging over disorder is equivalent to an averaging over the 
volume of a macroscopic system, an operation in which 
many microscopic details of the statics and dynamics be- 
come blurred. This might confer a mean-fleld character 
to the theory, where the contribution of the matrix would 
actually be taken into account at the level of a diffuse ef- 
fective localizing potential, with a possible loss of impor- 
tant local constraints. Unfortunately, this is a necessary 
step in order to develop a theory which is comparable 
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in complexity with the one for bulk systems, since it al- 
lows one to consider the system as homogeneous. Some 
progress has recently been made on the application of the 
MCT to inhomogeneous situations ^68*1, so it should be 
possible to relax the condition of homogeneity in order to 
study the dynamics of fluids confined in pores of simple 
geometry which are often preferred in simulation works, 
but there is no doubt that the more complex wavevector 
dependence of the resulting theory will make it harder to 
obtain numerical solutions of the equations. 

Moreover, beyond this purely technical aspect, the 
present formulation in terms of disorder averaged quan- 
tities has a practical interest as well. Indeed, many real 
porous media are disordered and most experimental tech- 
niques measure quantities which are averaged over the 
volume of a macroscopic sample and thus equivalent to 
the disorder averages considered by the MCT. So, the 
current theoretical setup seems well suited for direct com- 
parisons with experimental results. For molecular dy- 
namics simulations, however, since rather small systems 
are usually considered, it might be necessary to explicitly 
perform the disorder average over a representative sam- 
ple of matrix configurations before a comparison with the 
theory can be done. 

Altogether, in spite of the above merely technical is- 
sues, we believe that the present mode-coupling theory 
represents a valuable step towards a better understand- 
ing of the slow dynamics of confined glassforming liquids. 
Indeed, it remains rather simple and, since it is a mi- 
croscopic approach, it allows one to study in detail the 
effect on the dynamics of changes in the different ingredi- 
ents of a fluid-matrix model (fluid-fluid and fluid-matrix 
interactions, structure of the matrix). Thus, this pro- 
vides us with a tool to efficiently and thoroughly explore 
the phenomenology of dynamics in confinement. This 
is clearly illustrated by our findings for two very sim- 
ple systems with pure hard core interactions, which al- 
ready display new and nontrivial glass transition scenar- 
ios. Then remains the question of the validation of the 
theoretical predictions. Because the model of the QA 
mixture is quite simple and the theory makes detailed 
predictions, molecular dynamics studies should be able 
to give clear-cut answers. The presently available results 
look rather encouraging, but more simulation work is def- 
initely needed. 



APPENDIX A: REPLICA ORNSTEIN-ZERNIKE 
EQUATIONS 

For reference, we quote the replica OZ equations relat- 
ing the vario us p air correlation functions in QA binary 
mixtures [H, [3 [13, EBl ■ They read, in Fourier space, 



(Ala) 

+ n/c^ft-", (Alb) 
= + n^cl^h{^ + nfc{fh{f - n}c\h\, (Ale) 
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with cf^ = c'g + and h^^f ^ h'^ + h'^^. As usual, h 
and c denote total and direct correlation functions, re- 
spectively. / denotes the Fourier transform of / and the 
superscripts have the same meaning as for the structure 
factors (see Sec. II). 
Using the relations 
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(remember that S^^ = S'q + S",, hence = 1 
the OZ equations can be formally solved for the structure 
factors, leading to 
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APPENDIX B: DERIVATION OF THE 
MODE-COUPLING MEMORY KERNEL 
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In this Appendix, the derivation of the mode-coupling 
part of the memory kernel is outlined. This calculation 
is a direct extension of the one for bulk systems which is 
described in its most minute details in Ref. [l^ . 

The memory function in Eq. (jlip is defined as 



NfrnksT 



(Bl) 



where Rq{t) = exp[iQiCQit]iQiCgl^ is the projected 
random force obtained from the longitudinal momentum 



16 



density fluctuation g^{t) [5^. £ is the Liouville operator 
of the system and Qi is the complementary operator of 
the projector Vi which projects any dynamical variable 
onto the subspace spanned by Sp^ and g^. 

The calculation of the mode-coupling part of the ker- 



nel amounts to replacing i?q in Eq. (jBl|) by its projection 



V2Rq onto the subspace spanned by -Bq,k = 5pl^5p\ 
C. 



(1) 

q,k 



q-k' 



(5p{pq_k, and C^^^ = S p{^{p{^_-^) , where we intro- 



duce the projection operator V2 such that 



k 



'(i?qS-q,-k)Gq,ki?q.k + E E (^'l^-q.-k><f 



k 1,1' 



(B2) 



In this expression, we have anticipated that, within the 
mode-coupling approximation, V2 is diagonal in k and 
the subspaces spanned by the Bs and Cs are orthogo- 
nal. The prime in the first sum indicates that, to avoid 
double-counting, the wavevectors are assumed to be or- 
dered somehow and the sum is restricted to k < q — k. 

Gq_k and -ffq are normalization matrices insuring that 



V2Ba 



-2£>q,k = Bq,k and - k- 

In the course of this calculation, four-point density cor- 
relation functions are generated. Thus, in order to even- 
tually obtain closed dynamical equations, a factorization 
approximation is needed to express them as products of 
two-point density correlation functions. We follow the 
usual mode-coupling prescription and find that, within 
this approximation, the only non-vanishing four-point 
functions are given by 

^^/^.'^|q-k|0feW0|q-k|W, 

(B3a) 
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Specializing to t — and computing 7^2 ^q.k and 
^2^q k' results that 
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theorem, one finds 



qk 



2 "^Iq-kl 



q(q-k) 



/"^•^Iq-kl- 



(*^ffqCi'^._k) = ^qkBTN 
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Then, with the help of the extension of the convolution 
approximation [69'] to QA systems. 



(SpiSpi^Spi^^^) 



(B6a) 



(<5p^5/kP-q+k> = ^WN7nS',S^,S{^_^^, (B6b) 
(<5p^^/k(/q+k» = ^^/^,'^fe^|q-k| , (B6c) 



it comes 



(B7a) 



^ViCg^d'l_^) - tqkBT^7^,S'^S(^_^^ , (B7b) 



^riCg^d^l_^) = iqkBTNfSlS^^^y 



(B7c) 



The desired results are obtained by substracting the 
matching equations. 

We might now complete the explicit calculation of 
■P2-Rq- Indeed, if we define 



I'q.k = (-RqB-q,-k)Gq,k 
„(1) 



.(2) 



(21) 



(B8) 

%,ic = l^q^-q,-k;«q,k + l^q^-q.-k;^q,k ' (B9) 
<']c = (^qGi'i_k><k^ + (^qGil_k>i?S^ (BIO) 

such that 



It remains to express (i?qi3_q _k) = (^Qi'C^qB-q^-k) 

First, using Yvon's 



and (i?qGi'^ _k> = (^Qi-Cg^Gi 



^2i?q - E "'l^kSq^k + E <kC + <'kt^S, (Bll) 
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the above results immediately lead to 



q,k 

(2) 
q,k 
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q-k 1 q-(q-k) 1 



q Oj, 
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Note the vanishing of the fluid-matrix contribution. 

Eventually, injecting the resulting expression of 7^2 
into 



(^2-Rqft)^2-R-q) 

NfmkBT 



(B15) 



one obtains Eqs. (|13p and after a few elementary 
steps. 

Exactly the same calculation can be done with subsets 
of the above three quadratic variables. One then easily 
demonstrates that, if only Sp{^Sp^_^^ and '5pj^(pq_j.) are 
considered, the same equations are obtained, while, if 
one works with 5p{^Sp(^_i^ and Sp{^p^_j^, the equations 
of Ref. 371 are reproduced, which reduce to those of 
Refs. [29 JSO, [Ml in the zero fluid density limit. 



functional. Thus, it is especially suitable for our purpose, 
since it displays nonlinearities of the same degree as in 
the MCT. 

In the presence of a random magnetic field, just like 
QA mixtures are characterized by static frozen density 
fluctuations, the present spin-glass model is character- 
ized by static frozen local magnetizations (s;) which van- 
ish when the average over disorder (both on couplings 
and fields) is performed. Accordingly, we might define 
Ssi{t) = Si{t) — (si) and introduce three correlation func- 
tions, which are the analogues of S^, S^, and 4iq{t). They 
are, respectively, the two overlap functions 



1 ^ 



and 



1 ^ — 
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which obey Qb + Qc = Ij BLiid the normalized connected 
spin correlation function. 
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APPENDIX C: MEAN-FIELD SPIN-GLASS IN A 
RANDOM MAGNETIC FIELD 

A fruitful source of new theoretical developments on 
the physics of glassy systems has been the finding that 
the MCT provides an exact description of the equilib- 
rium dynamics of a certain class of rnean-field spin-glass 
models with multispin interactions [TOi, [7l|, [tI] . In this 
Appendix, based on Ref. [t^I (the interested reader is 
referred to this paper for details), we show that the 
equations describing the dynamics of such a mean-field 
spin-glass model in a random magnetic field reproduce 
the structure of those of the MCT for a fluid in a ran- 
dom environment. Thus, the correspondence between the 
two approaches still holds in the presence of an external 
source of disorder. 

We consider the fully connected mean-field spherical 
spin-glass model with three spin interactions in a random 
magnetic field, whose Hamiltonian is 



Hj\s] = - 



(1). 



(Cl) 



l<i<i<fe<JV 



l<i<A' 



where the spins Si are N real variables subject to the con- 
straint X)t=i ^1 — s-iid the random coupling constants 
Ji% fields j'f''' are uncorrelated zero mean Gaussian 
variables with variances 3J^/N^ and Ji/2, respectively. 
In zero field (Ji = 0), this model is known to generate 
the quadratic memory kernel introduced in Refs. [16j and 
[l3| as a one wavevector approximation to the bulk MCT 



The exact solution of the Langevin dynamics of the 
model using standard methods 73] provides us with an 
equation for the time evolution of (j){t) which reads, for a 
high enough temperature T, 



1 



1 



4>it) + —(f>{t) + — I dTm{t - T)(j){T) = 0, (C5a) 
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with 
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qlm^ + q. 
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A direct comparison of these equations with those of 
Sec. II immediately shows that both sets of equations 
exactly have the same formal structure. In particular, 
beyond the simple fact that the memory functions are 
similar polynomials of the dynamical correlation func- 
tions, the same parametrization in terms of connected 
and disconnected static correlations is found. 

We conclude with two short remarks. First, as for 
QA mixtures, we note that the external random field 
does not enter explicitly in the above equations (no Ji is 
present) . Because the model is exactly soluble, it is easy 
to understand how this happens. In fact, the details of 
the random field are only needed for the computation of 
the static correlations, through the equality 
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Once this calculation is done, they can be forgotten. The 
same is probably true for fluids in random environments 
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and this supports the assumption that the MCT equa- 
tions for QA systems actually have a wider domain of ap- 
plicability, provided they are expressed in terms of fluid 
quantities only. 

Second, it can be shown that Eq. (|C5P also describes 
the residual relaxation of the present spin-glass model in 
zero field, when the system is equilibrated in its low tem- 
perature phase and the total spin correlation function has 
a nonvanishing infinite time limit q [74]. Then, qt and 



qc simply have to be replaced by q and I — q. Thus, in 
the present model, the diff'erence discussed in Sec. Ill be- 
tween self-induced glassiness and the effect of a quenched 
random field does not seem to exist. This could have been 
anticipated, since by construction, fully connected mod- 
els are unable to capture phenomena which only show up 
in the wavevector dependence of the coupling constants 
in fluid systems. 
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